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We show that certain singularities of the Hamiltonian in the complex wave vector space can 
be used to identify topological quantum phase transitions for ID chiral topological superconduc¬ 
tors/superfluids in the BDI class. These singularities fall into the category of the so-called ex¬ 
ceptional points (EP's) studied in the context of non-Hermitian Hamiltonians describing open 
quantum systems. We also propose a generic formula in terms of the properties of the EP’s to 
quantify the exact number of Majorana zero modes in a particular chiral topological superconduct¬ 
ing phase, given the values of the parameters appearing in the Hamiltonian. This formula serves 
as an alternative to the familiar integer (Z) winding number invariant characterizing topological 
superconductor/superfluid phases in the chiral BDI class. 
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I. INTRODUCTION 


Exceptional points (EP's) are singular points in the parameter space of an operator at which two or more of its 
eigenvalues and eigenvectors coalesce m- If the operator is the Hamiltonian itself, EP ’s describe collapse of two 
or more energy eigenvalues at certain points in the parameter space similar to a degeneracy point, but with the 
important difference that the energy eigenvectors are not orthogonal to each other. The degeneracy of the eigenvalues 
with concomitant degeneracy of the eigenvectors gives rise to a whole host of non-trivial phenomena. For this reason, 
EP's have recently attracted enormous interest in the literature [201. Although the concept of EP’s is known in 
mathematics for many years, their application in physics has been mostly limited to open quantum systems with 
dissipation, appropriately described by non-Hermitian Hamiltonians im 

Topological superconductors [13] are systems characterized by a bulk superconducting gap, and, topologically pro¬ 
tected zero energy edge states known as MBSs (Majorana Bound States) described by second quantized operators 
satisfying the operator relation 7^ = 7. In the context of condensed matter physics, aside from being fascinating 
emergent non-elementary particles (which can be identified with their own anti-particles), MBSs obey Ising type 
non-Abelian braiding statistics mm potentially useful in implementing a fault-tolerant topological quantum com¬ 
puter mm- While MBSs have not yet been conclusively found in nature, they have been theoretically shown to 
exist in low dimensional spinless p-wave superconducting systems mi he], as well as other systems involving vari¬ 
ous heterostructures with proximity-induced superconductivity which are topologically similar to them |19H24| . In 
particular, the spin-orbit coupled semiconductor-superconductor heterostructure scheme has motivated tremendous 
experimental efforts with a number of recent works claiming to have observed experimental signatures of MBSs in 
zero bias tunneling experiments [25H2E] . More recently, experiments on ferromagnetic Fe-atom chains embedded on 
Pb superconductor substrate, have also seen tantalizing evidence of MBSs in spatially resolved scanning tunneling 
microscopy measurements |29| . 

Recent theoretical work [30, 31 has established that the quadratic Hamiltonians for gapped topological insulators 
and topological superconductors can be classified into ten topological symmetry classes, each of which is characterized 
by a topological invariant. The symmetry classification is important as it provides an understanding of the effects of 
various perturbations on the stability of the protected surface modes, such as MBSs. The (strictly ID) semiconductor- 
superconductor nanowire structure |22H24| , as well as the system of ferromagnetic atomic chains or nanowires deposited 
on Pb superconductor [22 [32 SSL are i n the topological class BDI. They are also known as chiral topological 
superconductors, described by an integer (Z) winding number topological invariant [34] that counts the number of 
protected zero energy Majorana modes at the individual edges. While the ID semiconductor-superconductor nanowire 
structure is topologically isomorphic to the ID spinless p -wave superconductor or Kitaev model |15| (which is in the 
chiral BDI class in the absence of symmetry breaking perturbations [34H36I ). the system of ferromagnetic atomic chain 
or nanowire embedded on Pb superconductor is isomorphic to the doubled (or time reversal symmetric) Kitaev model 
E2, which is also in the chiral BDI class with a Z invariant. In the presence of chiral symmetry breaking terms (say, 
for example, stray magnetic fields and/or magnetic impurities), the symmetry classification of these systems reduces 
to class D, described by a Z2 invariant, and the number of protected MBSs at any given end reduces to zero or one. 

The topological phases and quantum phase transitions within the BDI class topological superconductors are usually 
described in terms of the closing and re-opening of the single particle energy gap, and a winding number integer 
topological invariant associated with the bulk Hamiltonian with periodic boundary conditions that counts the number 
of protected zero energy end states localized at any given edge [32132132- In this paper, we describe the topological 
quantum phase transitions in the BDI class, and propose a generic formula to count the exact number of Majorana 
zero modes in ID chiral topological superconductors/superfluids based on the notion of EP' s. This quantity can 
exactly point out which topological phase we are considering, depending on the values of the parameters. In a recent 
work [22, the notion of the £P’s was discussed in the context of Majorana zero modes for the ID Kitaev model [Hi- 

Corresponding to a physical ID 2N x 2N Hamiltonian for a chiral topological superconductor, one can construct 
a non-Hermitian matrix by complexifying the momentum fc. In the chiral basis, the Hamiltonian has two N x N 
off-diagonal blocks, which is a property of a chiral topological superconducting system. The values of the complex 
fc, where the determinant of any one of the two off-diagonal elements vanish, are points in the complex k space 
where two (or more) eigenvalues of the complexified Hamiltonian vanish, and will be an example of the EP' s. These 
singularities are those special points where two or more repelling levels are connected by a square root branch point 
in the complex fc-plane [5|. More details can be found in Appendix [A] We will see that one or more EP 's collapse 
at the points of a topological quantum phase transitions in the parameter space. In other words, the determinants 
of both the off-diagonal blocks reduce to zero at the phase transition points. In the complex fc-plane, if one expands 
the off-diagonal Hamiltonian around a solution fc = kEP for an EP , the zeroth order piece is non-diagonalizable, 
as it consists of the determinant of one of the off-diagonal elements going to zero. The leading order correction is 
off-diagonal with determinant of both the blocks non-zero, and the entire matrix with these two non-zero off-diagonal 
blocks being diagonalizable. At a topological phase transition point, the the determinants of both the blocks of the 
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zeroth order matrix vanishes, and we are left with a completely diagonalizable Hamiltonian on expanding about that 
value of Rep- Elucidating the nature of a topological quantum phase transition in terms of EP’ s, we propose a generic 
formula in terms of the properties of the EP ’s to quantify the exact number of Majorana zero modes in a given chiral 
topological superconducting phase, given the values of the parameters appearing in the Hamiltonian. This formula 
serves as an alternative to the familiar integer (Z) winding number invariant [50} 1311 134] characterizing topological 
superconductor/superfluid phases in the chiral BDI class. 

T he p aper is organized as follows: In Sec.|TTJ we consider the Kitaev model [15] of ID spinless p-wave superconductor. 
Sec. Ill is devoted to the study of the quantum Ising chain with longer-ranged interactions [551140] . which, by a 
Jordan-Wigner transformation, maps on to the ID Kitaev model with longer range hopping and superconducting 
pair potential. We study a third model of the Majorana fermions in chiral topological ferromagnetic nanowires, 
with proximity-induced superconductivity [32} [33], i n Sec. IV All the three systems, for appropriate values of the 


parameters, can support chiral Majorana bound states at any given end. In Sec. [Vj we propose a generic formula to 
count the total number of zero modes in a particular chiral topological phase, given the values of the parameters of 
the Hamiltonian. We also provide the mathematical proof of why this quantity is related to the different topological 
phases. Lastly, we finish with a summary and outlook in Sec. m In Appendix [A] we review the definition of EP’ s. 


II. MODEL 1: KITAEV CHAIN 


Kitaev’s model of ID p-wave superconducting quantum wire can support Majorana zero modes at the ends m , 
depending on the value of the chemical potential fx. In the Bogoliubov-de Gennes (BdG) basis, the Hamiltonian is 
given by: 


Hi(k ) = (cos (k) — /x) cr y — A sin(fc) a x , 


with eigenvalues 


E = ± 




cos (k) — n) 2 + A 2 sin 2 (fc). 


(1) 

( 2 ) 


Rotating the basis, we write the Hamiltonian in the following off-diagonal form: 

A(k) = (n — cos (k)) + i A sin(fc), B(k ) = (/t — cos(k)) — i A sin(fc). (3) 

In the complex fc-plane, the EP’ s are given by 


( B ?t) T)' 


I u± 
r '-’AK 


= —i In 


—fx ± \/A 2 + /i 2 - 1 


A- 1 


} 


and kg K = —i In 


' M 


± \J A 2 + fx 2 — 1 


}■ 


(4) 


corresponding to A(k± K ) = 0 and B(kg K ) = 0 respectively. This is in conformity with our definition of EP ’s given 
in the introduction. Here the subscript “K” to the wave vector index k stands for Kitaev chain. 

Expanding around the EP’ s, we get: 

(k - k^ K ) , 

(k ~ 

(5) 


tfr(fc)- 

tfi(fc)- 


0 0 

B ( k AK ) 0 

0 A(kg K ) 

0 0 


0 

B'{k± K 

0 

b'(4 


A'{ki K 

0 

0 


The phase transitions occur at /z = ±1. For \fx\ < 1, we get a topological phase with one Majorana zero mode at 
each end of the chain. Moving from left to right along the //-axis, the transition are (1) from 0 to 1 chiral majorana 
mode as one crosses fx = — 1, and (2) from 1 to 0 as one crosses fx = 1. We note the following: 

1. At fx = —1, H\{k\ K ) = 0 and H\(k~^ K ) ^ 0, while H\{k^ K ) = 0 and Hi(kg K ) ^ 0. 

2. At fx = 1, H\(k\ K ) ^ 0 and H\{k~^ K ) = 0, while Hi(kg K ) ^ 0 and Hi{kg K ) = 0. 

3. At all other values of fx, H\{k^ K ) and {kj n< ) are non-zero. 
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FIG. 1. (Color online) (a) Blue region corresponds to Sr ( k AK ) < 0 and SI (k BK ) > 0. Purple region corresponds to SI {k AK ) > 0 
and SI < 0. (b) Blue region corresponds to SI (k AK ) > 0 and SI ( k BK ) < 0. Purple region corresponds to SI (k~^ K ) < 0 

and SI ( kg K ) > 0. 
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P 

FIG. 2. (Color online) f A ’ B (g, A) reproduces the topological phases of the Hamiltonian in Eq. |l]) in the (/r, A)-plane. 


So, one of the two EP 's collapses at the phase transition points, irrespective of whether we consider A{k) = 0 or 
B(k) = 0. _ 

show the dependence of the signs of SI (k AK bk) as functions of (/u. A). At g = —1, we find 


Figs. 


1(a) 


and 


1 (b) 


■c A k rsic(M, change signs. For instance, SI (k\ K (g, A)) < 0 for g < —1, and SI {k\ K (g, A)) > 0 for 

g > —1. This means that SI (k AK (g, A)) = 0 at g = —1, which in turn implies that E{k) = 0 has a solution for a 
real value of k (as E 2 (k ) = A(k) B(k)). Since E{k) = 0 for a real k signals a topological quantum phase transition 
(gap-closing) in this system, we find that in the PP-description, the topological phase transition is marked by the 
imaginary component of one of the FP’s going through zero. We find similar behaviour at // = 1 . These observations 
help us define the functions: 

f A ’ B (p A) = i sgn{Q (k\ KBK (g, A )) } ~ sgn{% [k\ KBK Oo, A)) } 


sgn{ SI (k AKBK (g, A)) } - sgn{ SI (k AKBK (g 0 , A)) } 


( 6 ) 


which capture the number of chiral Majorana zero modes in a given phase. Here, go is any value of g where we have 
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a non-topological or zero Majorana mode phase, i.e. g o ^ (—1,1). Fig. [2]shows the contourplot for f A ' B (g, A) in the 
{fi, A)-plane. In practice, one can work with either f A or f B to identify the topological phase. 

These EP’s were studied in an earlier work m in the context of topological phases for the Kitaev chain. The 


sgn{% [k\ K {g, A)) } + sgn{% (k AK (g, A)) } 


gives the number of Majorana 


authors proposed that the function ^ 

zero mode(s) in the various phases. However, we find that this definition works only for the 0 —» 1 Majorana 
fermion topological phase transition in the Kitaev model, and furthermore, when one chooses to work with the EP’s 
corresponding to A{k) = 0 (and not for B(k ) = 0). 


III. MODEL 2: ISING CHAIN WITH LONGER-RANGED INTERACTIONS 


The fermionized version of the ID transverse field Ising model, generalized to include longer-ranged spin-spin 
interactions, can support 0, 1 or 2 Majorana mode(s) at each end [38] • The corresponding BclG Hamiltonian is given 
by 


H 2 (k) = ( 
£(*0 = 2(1 


£(*0 *A 7 (fc)\ 

-iA T (k) -t(k) ) ' 

— Ai cos (A;) — A 2 cos(2 k)) , 


A j(k) = 2 (Ai sin(fc) + A 2 sin(2fc)) , 


(7) 


where Ai denotes the magnitudes of the nearest neighbour hopping and superconducting gap, and A 2 denotes the 
amplitudes of the next nearest neighbour hopping and superconducting gap. It has been shown [SB] that the presence 
of A 2 gives rise to two Majorana zero modes coexisting at the same end of the ID chain, for certain values of the 
parameters Ai and A 2 . In this system, multiple zero energy Majorana modes can coexist (and do not mix and split 
to finite energies), because the Hamiltonian is in the chiral BDI class with an integer invariant [31] .The eigenvalues 
of the Hamiltonian in Eq. ([TJ) are given by 


E I (k) = ±^e(k) + A 2 I (k). 


( 8 ) 


Since the Hamiltonian is chiral ( {H 2 (k), a x } = 0), by rotating the basis, we can rewrite it in the following off-diagonal 
form: 


H 2 (k) = 


0 Aj(ib) 
H/(fc) 0 


A I (k) = £(k)+iA I (k), B I (k)=£(k)-iA I (k). 


In the complex /c-plane, the EP’s are given by 


k AI — % In 


—Ai ± ^Xf + 4 A 2 


2 A 2 


} 


and k BI = —jin 


—Ai ± \/X\ + 4 A 2 


}■ 


(9) 


( 10 ) 


corresponding to Ai(k AI ) = 0 and Bi(k BI ) = 0 respectively. 

The phase diagram for this model has been derived previously [38]. Here we describe the topological phase diagram 
using the notion of the EP’s. Following the trajectories of the EP’s (derived by setting either Aj(k) = 0 or Bj(k) = 
0) as functions of (Ai,A 2 ), we find that along the phase transition lines from the 0 to 2 Majorana mode phase, 
H 2 (k AI bi) = 0, and hence both the EP’s collapse. However, along the transition lines from 0 to 1, or 1 to 2 
Majorana zero modes, only one of the EP’s collapses, i.e. either H 2 (k\ I BI ) = 0 or H 2 (k AI BI ) = 0 (but not both). 

show the dependence of the signs of 3 ( k AI BI ) as functions of (Ai, A 2 ). Let n denote the number 


Figs. 


3(a) 


and 


3(b) 


of Majorana zero mode(s) at each end. Along the n = 0 —>■ n = 2 phase transition lines, both 3 (k AI BI (Xi,X 2 )^j and 
k AT BI {Xi,X 2 )^ reduce to zero and go through a sign change. In Fig. j^J we show the contourplot for the functions 

f A ' B { Ai, A 2 ) = \ sgnl'A {k\ I BI { Ai, A 2 )) } - sgn{% [k\ I BI { A?, A°)) } 

+ sgn{% (k A i tB i( Ai, A 2 )) } - sgn{% {k Ai m { A?, X° 2 )^ } , 


( 11 ) 


where (A®, A°) is any point in phase space where we have a non-topological phase with no Majorana zero mode. This 
reproduces the topological phase diagram for the Ising chain with longer-ranged interactions found earlier 1 38] . 
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FIG. 3. (Color online) (a) Blue region corresponds to SI (k AI ) > 0 and SI < 0. Purple region corresponds to SI (fc+ 7 ) < 0 

and SI ( kgj ) > 0. (b) Blue region corresponds to SI ( k AI ) > 0 and SI (k BI ) < 0. Purple regions corresponds to SI (k AI ) < 0 
and SI ( kgj ) > 0. 



FIG. 4. (Color online) f A ’ B ( Ai,A 2 ), defined in Eq. (Ill, reproduces the topological phase diagram of the Hamiltonian in 
Eq. 0, found earlier [33 using the solutions of the BdG equations and the winding number topological invariant. Here, n 
indicates the number of protected zero energy MBSs at any given end. 


IV. MODEL 3: CHIRAL TOPOLOGICAL FERROMAGNETIC NANOWIRES 


We consider the strictly ID version of the Hamiltonian for a ferromagnetic nanowire embedded on Pb supercon¬ 
ductor [33] with a single spatial channel (i.e. no transverse hopping). In the momentum space, the BdG Hamiltonian 
becomes H = #3 (k) ’Pfc, where 


H 3 (k) = £ n (k) cr 0 T z + A s cr 0 + A p sin(fe) d ■ a t x +V-ctt 0 , £„(fc) = -2tcos(k) - fi 


( 12 ) 


Here k = k x is the ID crystal momentum, df/. = (cfc-j-, c/ 4 , cL fc ^, — c^_ k ^.) T is the four-component Nambu spinor which 
acts on the particle-hole (r) and spin (a) spaces, and V is the Zeeman field which can be induced by ferromagnetism. 
Also, A s and A p are proximity-induced s-wave and p-wave superconducting pairing potentials respectively, with d 
determining the relative magnitudes of the components of the p -wave superconducting order parameter A a p (a, f3 =t 
,4_). In our calculations, we use d = (1,0,0) and V = (0, V, 0). Furthermore, we set A s = 0, thus only considering 
p -wave pairing, which does not change the chiral BDI class of the Hamiltonian. 

The eigenvalues of the Hamiltonian are given by: 

Ei(k) =±y/(V + + A % sin 2 (k ), E 2 {k) = ±\J (V - £n{k )) 2 + A^ sin 2 (fc). 


(13) 
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We change the basis to transform the Hamiltonian in Eq. to the form: 


H (h\ — ( 0 hA ^ 

H ^ k) - l h B {k) 0 


where 

h A (k) = 


-£n(fe) +i Ap sin (k) 


V 


-v ’ '“><« = 


-£n(k) - iAp sin(fe) 

V 


- V 

€n(k) +iA p sin(fc) 


(14) 


This corresponds to the chiral basis where the chiral operator S = a x T y is rotated to the diagonal form diag(—1, —1,1,1). 
The upper block h A (k) and the lower block h B (k) of H 3 (k) give eigenfunctions of opposite chirality with respect to 
S. 

In the complex /c-plane, the EP 's for chirality solutions are given by 


k , 4,1 — 1 


i In | - 


V-fi± x /(V-^ + Al-AP 
2 1 + A„ 


} 


k A , 2 — 1 


i In | 


-V-n± y (V + Ai) z + A 2 - 4t 2 

2 1 Ap 


corresponding to det [ftx(fc)] = 0. Similarly, the EP 's for negative chirality solutions are given by 


kg 1 = —i In 


V- t i±\/(V- t i) 2 + Aj-4t* 


2 1 - A r: 


i kg 2 = —i\n 


-V — n ± yj(v + /z) 2 + A 2 - 4t 2 


2t — An 


}, (15) 


}, (16) 


corresponding to det [h B (k)] = 0. In either case, H 3 (k) becomes non-diagonalizable. 


Ei, 2 (k) 



FIG. 5. (Color online) Energy bands £4,2 (fc), given in Eq. (131, have been plotted in red and blue respectively as functions of 
p/t. We have used the values V = 1.5 1 and A p = 0.1 1. 


Choosing V = 1.5t,A p = 0.1 1, the energy bands are shown in Fig. ^ as functions of /z/i. So, there are phase 
transitions at /i = {—3.5 f, 0.5 f} for h 2 {k), and at fi = {—0.5 1, 3.51} forTii(fc). Moving from left to right along the 
/z-axis, the transitions are (1) from 0 to 1 chiral Majorana mode as one crosses /i = — 3.5 f, and (2) from 1 to 2 chiral 
Majorana modes as one crosses /jl = —0.5 1. Let us examine the EP' s corresponding to det [ h A (k )] = 0: 

1. At /z = -0.5 1, det [h B (k] J^)] = 0, while det [h B (k A x )] and det [h B (*£ 2 )] are non-zero. Hence, H 3 (k\ -, ) is 
diagonalizable indicating the collapse of the EP at k = k A1 . 

2. At /z = 3.5 1, det[h B (k A1 )] = 0, while det [h B (k A1 )] and det \h B (k% 2 )] are non-zero. H 3 (k A1 ) is thus 
diagonalizable indicating the collapse of the EP at k = k A1 . 

3. At fi = —3.5 1, det [h B (k A2 )] = 0, while det [h B (k A2 )] an d det [h B (k% 1 )] are non-zero. H 3 (k A2 ) is thus 
diagonalizable indicating the collapse of the EP at k = k A2 - 

4. At n = 0.5£, det [h B (k A2 )\ = 0, while det [h B (k A2 )] and det [h B (k% 1 )] are non-zero. H 3 (k A2 ) is thus 
diagonalizable indicating the collapse of the EP at k = k A2 . 












sgn{lm{*J'i)} 


sgn{lm(*^)} 


■f 


■+ 


-5 -4 -3 -2 -1 12345 




-5 -4 -3 -2 -1 1 2 3 4 5 




(a) 


(b) 


(c) 


FIG. 6. (Color online) Parameters: V = 1.5 1, A p = O.lt. (a) Blue and dotted red lines correspond to s<?n{A (k\ -J j and 
sgn{ A ( k^ x ) } respectively, (b) Blue and dotted red lines correspond to sgn{£s (k\ 2 ) } and sgn{Q (fcjj 2 ) } respectively, (c) 
fsum(g) giving the count of the chiral Majorana zero modes as functions of g/t. 


Similar observations hold if we consider the EP's corresponding to det [hg{k)] = 0. So, one of the four EP 's collapses 
at each phase transition point, irrespective of whether we decide to examine the det [/^(A;)] = 0 or det [hg (k)] = 0 
solution 
Figs. 


6(a) 


and 


6(b) show the dependence of the signs of 3 ^ and 3 on g/t. Generalizing the functions 

defined in Hqs. |6]) and to the present case of four EP' s, we define: 

1 


f£Im (m) = 2 I [ S 3 n { 3 ( fc A/bA v)) } ~ S 9 n {% ( k A/B, r(W>)) } 


s=± r—1,2 


(17) 


which captures the number o f chira l Majorana zero modes in a given phase. Here, go is any value of g where we have 
a non-topological phase. Fig. 6(c) shows the plot for fsum (m)- 


V. GENERIC FORMULA 


A generic ID 2 N x 2 N chiral Hamiltonian H chirlij {k) can support multiple Majorana zero modes at each end of the 
fermionic chain. Let k = k 3 EP ( for j = 1,2,..., n ep ) be the solutions of the EP 's corresponding to the vanishing of 
the determinant of any one off-diagonal block, after H ( ±j ra i has been rotated into the off-diagonal form and k has been 
promoted to a complex number. From the study of our models, we observe the following: 


1 . 


In Sec. Ill sgn{kj {k\ K {{pi})) } = sgn{kj (k AK ({pi})^ } for a phase with n = 1 zero mode. Furthermore, 
sgn{ks {k^ K ({pi})) } and sgn{ks (k^ K ({pi})) } have opposite signs for the non-topological phases. 


2. In Sec. 


Ill 


sgn{3 (k~Xi{{Pi})) } = s gn{ks (k^j({pi})) } for a phase with n = 0,2 zero modes, and these 
signs are opposite to each other for the n = 0 and n = 2 cases. Furthermore, sgn{^s (k\j({pi})) } and 
sgn{ks (k^ r ({Pi})) } have opposite signs for a phase with n = 1 zero mode. 


3. 


In Sec. |IV[ one of the four EP 's undergo sign change at each phase transition point, and the sign dependence 
can be examined from Figs. |6(a~)1 and [6(b)] 


Here we discussed only the results for the upper off-diagonal block. But similar results hold for the lower one. 
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Hence, we propose the following generic formula: 

_| flep 

/(M) = 2 | H ( fc i;p(fe})) } - sgn{$s (k 3 EP {{p?}y) } 


i=i 


(18) 


where (for j = 1,2,... n ep ) are the DP solutions found by solving any one of the off-diagonal blocks, {pi} is the 
set of parameters appearing in the expressions for k EP , and {p] 1 } are their values at any point in the non-topological 
phase. 

We emphasize that we can apply this formula for any chiral system in ID, i.e. for any chiral topological super¬ 
conductor in the BDI class with an arbitrary integer topological invariant Z. One example is the the case of several 
nanowires coupled by a transverse hopping term. For the BDI class in ID, MBSs appearing at one end are of the same 
chirality, which at the same time characterize the winding number in that phase. Hence it is sufficient to consider the 
DP’s corresponding to the solutions of one of the off-diagonal blocks in order to apply Eq. (18). 

We now prove why one (or more) 3 (^k^ P J changes sign at a topological phase transition point, characterised by a 
set of parameters {pi = p\}. We consider a ID 2TV x 2 N Hamiltonian 


n (h\-( 0 MMpJ) 

U B ( fc;fe }) o 


(19) 


after rotating it to the off-diagonal form, such that the N x N matrices Ha and h E satisfy h B (k , {p.;}) = k' A {k, {pi}) 
for real (physical) values of k. Let 


det [h A (k,{pi})\ = a(k, {pi}) - ib(k,{pi }), 

where the parameters {p*} are real, and the functions (a, b ) are real for real k. Then it follows that 

det [h B {k,{pi})] = a(k,{pi}) +ib(k,{p z }). 


( 20 ) 


( 21 ) 


Let us consider the DP’s at k = kEP in the complex fc-plane corresponding to det [h A (k, {p*})] = 0 => a(kEP , {pi}) = 
ib(k E p,{pi})- Let 


k E p = x{{Pi}) + w({Pi}), 


( 22 ) 


where ( x,y ) are the real and imaginary parts of k E p as functions of {pi}. Then we must have 

a(k EPl { m }) = /r en (x, y) + if° a dd {x, y ), b(k EP , {pj) = / 6 eoen (®, y ) + ^f 0 b dd {x^ y) , (23) 

such that f™ b en {x,y) and fa d b (x,y) are even and odd functions of y respectively. Needless to add that they must be 

real functions too. Now the constraint a(k E p , {p,;}) = i b(k E p , {p;}) translates into the equations 

f e a ven (x,y) = f° b dd {x , y ), f° a dd {x , y) = /r"(z, y) , (24) 

whose solutions give (x,y) as functions of {pi}. At a point {pi = p(} in parameter space, where both a ± ib vanish 
(i.e. det [h A ] = det [ h B \ =0), we must have 

a{{pVi) = H{Pi}) =0 => fa dd (x,y)\{ Pi =pi} = f b dd (x,y)\ {pi=p t } =0 => y{{p\}) = 0 (25) 

as one of the solutions for p({p-}). This is clearly a topological phase transition point, as det [H 0 d(k)\ vanishes for a 

real value of k, which indicates one or more energy eigenvalues going to zero. At this point, the DP also collapses as 
the Hamiltonian becomes diagonalizable. For a point {piiP--^} close to {p\ ,p*_^ 1 }, we have the expansion 


2/({PuP^i}) 


\pi - Pi I dp\y{{Pi }) if Pi > p\ , 

-\pi - p\ I dpt y{{p\}) if Pi < p\ , 


(26) 


for the solution y({p\}) = 0. Hence, this solution for 3 (k EP ) undergoes a sign change on crossing the phase transition 
point. 

Our counting formula is supplementary to other counting methods studied in the literature, like winding number 
calculation [SHUSHES!, scattering matrix approach [41] 02] and gradient expansion [43]. The advantage of using our 
formula is that one has to solve only for the zeroes of the determinant of one off-diagonal block in terms of complex 
k. One need not perform any integral or other lengthy computations to find the number of MBSs characterizing the 
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topological phase. Hence, in spite of conveying the same physics, we believe this method is slightly easier and more 
convenient to implement than the other standard methods. 

From our entire analysis, we can intuitively understand that these EP solutions in the complex /c-plane can be 
mapped to the MBS wavefunctions in the real space with open boundary conditions. At any one end of a chiral ID 
system, only MBS wavefunctions of a definite chirality can appear, which is reflected by the fact that we need to solve 
for EP 's only for one off-diagonal block and this gives us the number of MBSs at each end of an open chain. Again, 
the fact that 3 (k EP ) changes sign at a topological phase transition is tied to the fact that an MBS mode appears or 
disappears on moving from one phase to the other. One can prove the correspondence exp (ikppx) -H- exp(— zx), 
where exp (—z x) represents an exponentially decaying MBS wavefunction at a distance x from the edge. Thus the 
imaginary part of momentum at the exceptional point indeed gives us the exponential decay of the corresponding 
MBS in real space. A rigorous proof of this bulk-boundary correspondence has been provided in subsequent works 

[3311351 . 

We note that we cannot apply this EP formalism to count the Majorana zero modes for ID Hamiltonians in class 
D. This is because a chiral symmetry operator does not exist and the Hamiltonian cannot be unitarily rotated into 
the block-diagonal form. It can be represented as: 


H D (k) = cos 6 


0 h A (k, {pi}) 

hp(k, {pi}) 0 


+ sin 9 


I 0 \ 

0 —I ) ’ 


(27) 


with a non-vanishing value of sind. In such a situation, the EP solutions for the Hamiltonian will not correspond 
to the vanishing energy eigenvalues in the complex fc-plane. For class Dill systems, the Majorana fermions at one 
end of the edge are not of the same chirality, but correspond to Majorana Kramers pairs (MKPs) which are doubly 
degenerate Majorana zero modes. So to count those modes, we have to apply some additional criterion to the counting 
formula. The counting formula for the cases of class D and Dill systems has been discussed in subsequent works 

mm- 


VI. CONCLUSION 

We have derived a generic formula for counting the number of Majorana zero modes for a ID chiral topological 
superconductor/superfluid. First we write the Hamiltonian in the basis where the chiral symmetry opertaor is diagonal. 
In this basis, the Hamiltonian consists of two off-diagonal blocks. The solutions, k EP (for j = 1,2, ...n ep ), for the 
exceptional points in the complex momentum space, can be obtained in terms of the parameters of the Hamiltonian 
by setting the determinant of any one off-diagonal block to zero. Our formula is based on the evolution of these EP's 
in the complex fc-plane as functions of the parameters. The count of the MBSs is encoded in the signs of 3 (k EP J ■ 
At an EP , by definition, the complexified Hamiltonian is non-diagonalizable, due to one of the eigenvectors having 
a vanishing norm. However, at a physical phase transition point, at least one of the k J EP s becomes real making the 
Hamiltonian diagonalizable. This is due to the vanishing of the determinants of both the off-diagonal blocks, thus 
signalling the collapse of the corresponding EP. 
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Appendix A: Review of Exceptional Points 

“Exceptional Points” are branch point singularities in the parameter space of a matrix, at which at least two or 
more eigenvalues coincide. These are different from the familiar degeneracy observed for a Hermitian operator because 
of the fact that both the eigenvalues and the corresponding eigenstates coalesce [5]. Since a Hermitian operator has a 
complete set of eigenstates, EP 's can occur only in the spectrum of a non-Hermitian operator. Hence, if we analytically 
continue a real parameter of a Hamiltonian to complex values, the resulting system becomes non-Hermitian and we 
may find values of that complexified parameter where EP’s appear. Let us demonstrate the appearance of the EP 's 
for the simplest case of a 2 x 2 matrix: 



(Al) 



12 


The eigenvalues of TL are given by: 



jfi + £2 + 7( £ i + <7>) ± 


£2 + £i7 - «27) 2 + 4 7 2 r7i?72 


}■ 


(A2) 


which coalesce at 


(£1 - + ei7^ e 2 7) 2 +47 2 77!77 2 = 0 => 7 = j ± = ^ = • ( A3 ) 

H e l ^ e 2j ± V 7 ?!^ 

At these points, TL is a non-diagonalizable non-Hermitian matrix, which is manifested by the fact that there is only 
one linearly independent eigenvector (instead of two) and with a vanishing norm. The bi-ortlrogonal system for a 
non-Hermitian operator gives the right and left eigenstates proportional to: 

, ( , - >1) for 7 = 7+; 

Vvi m 

, ( lT]2 ,1) for 7 = 7_ , (A4) 

Vvi m 

which clearly have zero norm. Generalization to a higher dimensional matrix is straightforward. 




